Ratchet effect in surface electromigration: 
smoothing surfaces by an ac field 
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We demonstrate that for surfaces that have a nonzero Schwoebel barrier the application of an ac 
field parallel to the surface induces a net electromigration current, that points in the descending step 
direction. The magnitude of the current is calculated analytically and compared with Monte Carlo 
simulations. Since a downhill current smoothes the surface, our results imply that the application 
of ac fields can aid the smoothing process during annealing and can slow or eliminate the Schwoebel 
barrier induced mound formation during growth. 

PACS numbers: 68.35.Ct, 68.35. Md 



Growing epitaxial films with smooth surfaces is one of 
the ongoing challenges of the thin film community. How- 
ever, this goal is hampered by a series of basic physical ef- 
fects that lead to the development of unavoidable surface 
roughness during growth. In particular, there is abun- 
dant experimental and theoretical evidence that during 
deposition the diffusion bias generated by the Schwoebel 
barrier (see Fig. |l|) results in a net uphill current, which in 
turn leads to the formation of mounds and to a fast and 
unwanted increase in the interface roughness jjj . As Fig. 
|l| demonstrates, the Schwoebel barrier introduces spa- 
tial asymmetry in the otherwise symmetric lattice poten- 
tial. Interestingly in the recent years has been recognized 
that in such periodic and spatially asymmetric systems 
(ratchets) non-equilibrium fluctuations can rectify Brow- 
nian motion and induce a nonzero net current |Q. This 
fluctuation driven transport is believed to play an essen- 
tial role in the operation of motor proteins or molecular 
motors, and might result in new separation techniques 
In this paper we propose the first nano-scale appli- 
cation of this transport mechanism based on the atomic 
electromigration on vicinal surfaces induced by alternat- 
ing electric fields. We demonstrate that the Schwoebel 
barrier induced asymmetry in the lattice potential can be 
used to generate a downhill current, aiding the smoothing 
of surfaces during growth or annealing. 

Atom diffusion on crystal surfaces is a thermally ac- 
tivated process: atoms can hop from their position to 
a neighboring one by overcoming a potential barrier 
AE. The hopping rate is given by the Arrhcnius law 
k = V{) cxp(— AE / ksT) , where T is the temperature and 
vq is the vibration frequency of the surface atoms. Fig. 
|l| illustrates the lattice potential of a vicinal surface that 
consists of long flat terraces separated by monatomic 
steps. The barrier height for diffusion on a flat surface 
is denoted by E . Near a step atoms form additional 
lateral bonds of energy E\ with the step atoms, lead- 
ing to a deeper potential valley. Finally, jumping over a 
step requires passing an additional potential barrier, the 
Schwoebel barrier, Eh M. 



For most metals and semiconductors the otherwise ran- 
dom surface diffusion of the atoms can be biased by an 
external electric field applied parallel to the surface, a 
phenomenon known as surface electromigration Q . The 
effective force, F, acting on the surface atoms is propor- 
tional to the field E, F = ZeE, where e (> 0) is the 
elementary charge and Z is the effective charge number 
which consists of two terms, Z = Z^ + Z w . 





FIG. 1. Schematic illustration of the cross section of a 
vicinal surface containing two monatomic steps (A) and the 
asymmetric lattice potential (B) experienced by an atom dif- 
fusing on this surface. Note that if the Schwoebel barrier 
Eb = 0, the lattice potential near a step is spatially symmet- 
ric, while this symmetry is broken for Eb 7^ 0. A downhill 
(uphill) current points to the right (left) on this figure. As 
(C) illustrates, if there is a mound on the surface at t — 0, a 
downhill current indicated by the arrows will tend to decrease 
its height, i.e. it smoothes the film by the reorganization of 
the material on the surface. 

The "direct" term Z^ (> 0) is associated with the elec- 
trostatic interaction between the atom and the electric 
field, while the "wind" term Z w (< 0) is generated by 
the scattering of the current carrying electrons on the 
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surface atoms. The competition between these two terms 
can result in either positive or negative effective charge 



The constant electric field induces a surface current 
parallel to the field. However, it is expected that a pe- 
riodic field with zero mean does not create a net current 
over a full period, since it simply amounts to biasing in 
an equal manner the motion of the atoms in two oppos- 
ing directions. In contrast, next we demonstrate that 
for systems that have a nonzero Schwoebel barrier the 
asymmetry induced by the barrier rectifies the diffusion 
process, generating a net current along the surface, even 
if the external field has a zero mean value. Most im- 
portant, the induced current is always downhill, i.e. it 
points towards the descending step direction, indepen- 
dent of the step orientation or the effective charge. Since 
the downhill current acts to smooth the surface, it has 
the potential to accelerate the smoothing process during 
annealing and to slow or eliminate the Schwoebel barrier 
induced mound formation during growth. Consequently, 
this nano-scale ratchet effect can have important techno- 
logical applications for thin film growth. 



Consider a vicinal semiconductor surface with terraces 
of equal width w in an alternating external electric field 
(with zero mean) perpendicular to the steps and paral- 
lel to the surface. The lattice potential of this system 
is periodic and spatially asymmetric, and the alternat- 
ing field acts as a zero-mean non-equilibrium fluctuating 
force. In the presence of such field (or ac current) a 
directed net flow of the diffusing atoms on the surface 
is expected from the theory of fluctuation driven trans- 
port (|]. A simple explanation of the effect responsible 
for this current is given in Fig. ^. To evaluate the net 
current we first consider the motion of a single atom on 
a fixed vicinal surface with terraces of width w. Since 
the diffusion parallel to the steps is not affected by the 
electric field we can neglect the transverse direction and 
consider the motion in one dimension, perpendicular to 
the steps. The one-dimensional lattice potential for the 
diffusing atom is similar to the potential of Fig. [ij but it 
is periodic with period w, each period containing £ val- 
leys and barriers, where £ — w/a is the dimensionlcss 
size of the terraces measured in the units of the lattice 
constant a. Associated with the steps, every (£ — l)th 
Eq barrier is followed by a higher Schwoebel barrier Et> 
and a deeper valley E\ representing binding to the step. 
The system can be reduced to one period of the potential 
with periodic boundary conditions, in which the motion 
of the atom is described by the master equation 
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FIG. 2. The effect of an uniform electric field on the po- 
tential landscape experienced by an atom diffusing on the 
surface. The line (A) corresponds to the lattice potential in 
the vicinity of a step for E = (see also Fig. Up). Assum- 
ing Z > 0, a large enough field to the right (E > 0) 'tilts' 
the potential landscape, eliminating the minima of the po- 
tential, biasing the adatom motion to the right (see (B)). On 
the other hand, as (C) illustrates, an inverse field (E < 0) 
of the same magnitude does not eliminate the minima in the 
vicinity of the step (marked by an arrow), locally trapping 
the particle. Consequently, the current to the left in (C) is 
smaller than the current to the right in (B), and a periodic 
external field will result in a net current to the right, i.e. a 
downhill current. For smaller fields there are minima in both 
cases, but the minimum corresponding to (C) is deeper than 
that of (B), again resulting in a more efficient trapping of the 
atoms moving uphill, and to a net downhill current. Note 
that the direction of the net current is determined only by 
the step orientation, and not by the sign of Z. 

(i) 

where 1 < i < £ with periodic boundary conditions, Pi(t) 
denotes the probability that the atom is situated in the 
«th valley at time t. k~ and kf denote the hopping 
rates from the ith valley to the left and to the right, 
respectively, where fcr(i) = u c -[ AE i + SE ^}/ k * T , and 

k+{t) = v oe ~l AE >~ SE ^y k » T , Ai?r - E , AE+ = 
E + E h , AEz =E + E!+E h , AE+ =E +Ei, and 
AEj = AE+ = E for 3 < j < £ are the heights of the 
barriers of the lattice potential. Finally, 5E(t) — F(t)a/2 
is the change in the barrier heights associated with the 
force F(t) — ZeE(t) induced by the electric field E(t). 

If the period of the electric field is much larger than 
the relaxation time r ss £ 2 ey^(Eo / k-oT) / of the prob- 
ability distribution Pi (t) , Pi (t) is "slaved" by the electric 
field and the solution of the master equation reduces to 
the stationary solution with constant electric field E. If 
the period of the field is much smaller than r, Pi(t) has 
no time to adjust to the field, resulting in a very small 
net current. Thus /„ = I/r is the upper limit of the 
frequency, since exceeding it the current goes to zero. 
Therefore, in the following we assume that the frequency 
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/ of the alternating field is below the upper limit /„, so 
that the stationary solution of the master equation ap- 
plies. In this case Eq. ([!]) reduces to a system of {£ — 1) 
independent linear equations with the normalisation con- 
dition 5Zi=i Pi = 1) gi ym S the £th equation. The solution 
provides the stationary current as 

J st (F) = k+(F)P 1 (F) - fc 2 -(F)P 2 (F) = 
= 2 smh(£SE/k B T) v e~ Eo/kBT x 



1 e 



(l-t)5E/k B T 



sinh(«£/fc B T) 

sinh((5 J B/fc B T) 



(2) 



For simplicity we restrict the calculation to the case 
when the field is as a symmetric square wave, i.e. the 
force alternates between +F and —F at constant inter- 
vals. In this case, when / <c /„, the net current Jo can 
be calculated from (|) as J = [J s t{+F) + J st (~F)}/2. 

If the interaction between the diffusing atoms is ne- 
glected, the average number of atoms that are able to 
move (i.e., that are not part of a terrace) is 1+N a in each 
period. The first term indicates the edge of the terrace 
and the second term, N a = (£— 1)/[1 + exp^i/feeT)], is 
the average number of diffusing atoms (adatoms) on the 
surface. Thus the net particle current is 



J = MF){1+N a ) = J a (F)[ 1 + 
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(3) 



Eq. (|^) provides the downhill current generated by the 
interplay between the ac field and the Schwoebel bar- 
rier. The terrace size dependence of this current for two 
different temperatures is shown in Fig. ||A. 

The previous calculation, while correctly describes the 
nature and the qualitative features of the net current, 
neglects the atom-atom interaction and the step fluctua- 
tions. Since the source of atoms are the steps (adatoms 
detach from step edges) , the step length is not fixed, but 
it fluctuates. To incorporate these effects and to check 
the validity of the analytical prediction, we performed 
Monte Carlo (MC) simulations with activated diffusion 
along the surface. In the simulations we start with a se- 
ries of steps of length £. Every surface atom that has less 
than two neighbors in the same layer is allowed to diffuse 
with a probability P ~ exp(— AE/ksT). In the presence 
of an electric field SE the activation energy is reduced or 
increased with SE, depending on whether the direction 
of the hop is along or against the field. In the simulations 
we use a periodic square wave with frequency / <C /«, 
and determine the net current defined as the number of 
excess hops in the downhill direction in the unit time, 
normalized to the system size. 

As Fig. ||A indicates, the net current obtained in the 
Monte Carlo simulation qualitatively agrees with the pre- 
diction (||) of the master equation. Furthermore, there 



is excellent quantitative agreement at low temperature 
(900 K) while (|) overestimates the current at higher T. 
Note that the fit does not require any fitting parameter. 

To pinpoint the temperature range for which the an- 
alytical solution provides a good approximation, in Fig. 
Hb we plot the temperature dependence of the net cur- 
rent for £ — 20 and £ — 200. We find excellent agreement 
between the theory and the MC simulations for T < T* 
(~ 1000 K), while for T > T* the current predicted by (§) 
systematically exceeds the numerical result. Indeed, in- 
creasing the temperature also increases the adatom den- 
sity on the surface, and consequently the role of the ne- 
glected atom-atom interactions also increases The 
value of T* depends only on the energy barriers Eq,E\ 
and Eh- In some materials the bonding energies might be 
higher than the quoted values, which further increases T* 
and the range of applicability of the analytical solution. 
Note that the agreement between the analytical predic- 
tions and the MC simulations is better for £ = 200 than 
I = 20, underlying the larger impact of step fluctuation 
on shorter steps [pj. 
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FIG. 3. (A) The ac field induced downhill current as a 
function if the average step size I. The lines correspond to 
the analytical results (^) and Q for T=900 K (continuous 
line) and T=1200 K (dashed line), while the circles (900 K) 
and triangles (1200 K) are the currents predicted by the MC 
simulations. (B) The temperature dependence of the net cur- 
rent for average step sizes I = 20 and i — 200. The continuous 
(I — 20) and dashed (£ = 200) lines correspond to the ana- 
lytical predictions, while the circles (I = 20) and triangles 
(I = 200) are the result of the MC simulations. In the sim- 
ulations we used the parameters Eq = 0.3 eV, E\ = 0.6 eV, 
E b = 0.15 eV, Z = 0.5, E = 10 s V/m, and v a = 10 13 s" 1 . 
The system size is L — 2000 and the results were averaged 
over 20 independent runs. 
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Since in the typical electromigration experiments || 
SE is much smaller than fcgT, we can expand the current 
Jo(-F) into Taylor series in terms of (SE/ksT), obtaining 
that the net current is a second order effect, proportional 
to (Z8E) 2 . The MC simulations confirm this prediction, 
providing a quantitative expression for tuning the current 
with E P|. 

In conclusion, we have demonstrated both analytically 
and numerically that the Schwoebel barrier, in the pres- 
ence of a periodic external electric field, leads to a down- 
hill current. Since most metal and semiconductor sur- 
faces do have a nonzero Schwoebel barrier and display 
electromigration, we expect that the appearance of such 
a net current is relevant for a large class of technologi- 
cally important materials. Thus the application of an ac 
current during either growth or annealing can lead to a 
nontrivial smoothing effect, and aid the growth of smooth 
surfaces. This consequence of the ratchet effect can thus 
have important practical applications in the growth and 
processing of high quality thin films. 
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